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ABSTRACT 

The  reaction  operator  formulation  for  the  exact  solution  of  the 
SchrOdinger  equation  is  used  with  a  truncated  basis  set  to  obtain 
approximate  solutions.  The  relationship  between  this  truncated  reaction 
operator  formalism  and  the  Rayleigh-Ritz  variational  method  is  emphasized 
and  shown  explicitly.  The  truncated  reaction  operator  and  its  matrix 
elements  are  discussed  in  general;  computed  and  discussed  for  a  simple 
example.,  the  helium  atom^  using  a  three  membered  basis  set.  It  is 
shown  that  the  reaction  operator  can  be  replaced  by  a  function  which  we 
call  the  "effective”  perturbation.  This  function  has  very  fundamental 
significance  and  may  lend  itself  to  accurate  empiricism.  Expressions 
are  given  comparing  the  "effective”  perturbation  with  the  perturbation 
and  approximations  to  this  function  for  the  helium  atom  are  computed. 
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TRUNCATED  REACTION  OPERATORS 


1°  Introduction 

1  2  3  4 

Watson  ^  Brueckner  and  Levinson  5  LBwdin  and  others  have  derived 

the  reaction  operator  formulation  for  the  exact  solution  of  the 

SchrUdinger  equation,,  This  formulation  can  be  used  with  a  truncated 

basis  set  to  obtain  approximate  solutions  of  the  SchrUdinger  equation.. 

One  of  the  purposes  of  our  present  discussion  is  to  emphasize  and  show 

explicitly jj  in  the  simplest  possible  manner^  the  relationship  between 

the  approximate  solutions  obtained  with  a  truncated  reaction  operator 

formalism  and  the  solutions  obtained  by  the  Rayleigh-Ritz  variational 

method*  The  approximate  reaction  operator  and  its  matrix  elements  are 

discussed  in  general  and  in  connection  with  a  simple  example^  i,e0  the 

helium  atom^  using  a  three  membered  basis  set* 

The  reaction  operator  can  be  replaced  by  a  function^  V  ^  >  which 

we  call  the  "effective”  perturbation.  It  appears  to  us  that  this 

function  has  very  fundamental  significance 0  Ultimately  after  we  have 

obtained  a  backlog  of  numerical  and  theoretical  experience  with  the 

determination  of  the  "effective"-  potentials  for  atomic  and  molecular 

problems 9  we  may  be  able  to  accurately  empiricize  V  ^  9  thus  obtaining 

good  values  for  the  energy  and  other  physical  observables  of  these 

systems . 

Firsts  consider  that  part  of  the  reaction  operator  formulation 
for  the  exact  solution  of  the  Schrftdinger  equation  pertinent  to  this 
paper.  We  fix  our  attention  on  the  effect  of  the  perturbation  V  on 
the  q-th  unperturbed  state.  The  q-th  state  is  considered  to  be  non~ 
degenerate  and  is  not  necessarily  the  ground  state.  The  Hamiltonian 
for  the  unperturbed  system  is  and 

^  K.  M.  Watson^  Phys.  Rev.  89^  575  (1953). 

2 

K.  A.  Brueckner  and  G.  A.  Levinson Phys.  Rev,  97,,  1344  (1955), 

3  P.  0.  LBwdin,  J.  Math.  Phys.  3,  969  (1962). 

4 

A  more  complete  list  of  references  is  given  in  K.  Kumar  ^  "Perturb a~ 
tion  Theory  and  the  Nuclear  Many  Body  Problem"  (Nor th-Hol land 
Publishing  Co.  9  Amsterdam^  1962). 
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Hofk  =  Ek°^k  >  k=l, 2,  •  •  •  ,  (1) 

where  and  Ek.°^  are  norma^:*-ze^  eigenfunction  and  eigenvalue, 

respectively,  for  the  k-th  unperturbed  state.  For  the  perturbed 
system,  the  Hamiltonian  H  ,  is  given  by 

H  -  H  +  V  (2) 

o  v  ' 

The  reaction  or  energy  shift  operator,  t  ,  is  defined  by 

\  ■  E$o)  +  <s9?,ltl f£>  <3) 

where  E^  is  the  exact  energy  for  the  perturbed  system.  Let  us  define 
two  projection  operators,  0  and  P  ,  such  that 

OX-?, 

and 

P  =  1  -  0 

where  X  is  a  trial  function.  Using  these  definitions  it  can  be 

shown  that  the  hermitian  operator  t  ,  associated  with  V  ,  H  , 

3  0 

and  the  state  under  consideration  q  ,  is  given  by 

t  -  V  +  VT  t  (6) 

o  ' 

where 

To  =  p  [<X-0  +  P(Eq  -  Hq)pJ  _1p  .  (7) 

Here  oc  is  an  arbitrary  number  which  insures  the  existence  of  the 
inverse  operator  in  Eq.  (7) . 

3 

The  exact  wave  function  for  the  perturbed  state,  i|r  ,  is  given  by 

*q  =  (8) 


(4) 

(5) 


3 

where  W  is  the  wave  operator  associated  with  the  q-th  unperturbed 
state; 


W  -  1.  +  T  t  =  (1/V)  t 


(9) 


Here 

<?,'*,>  -1 


(10) 


and  hence  the  wave  function,,  \|r  ^  has  not  been  normalized 

One  may  define  an  "effective"  perturbation^  > 

replaces  the  operator  t  in  the  above  formulation; 


to  unity, 
which 


W»  -  (t  ■  (11) 

Clearly  the  exact  energy  and  wave  function  for  the  perturbed  system 
may  be  written  as 


\  -  W’>l?q>  <l2> 

and 

+q  =  <ve£f  (q)/V)<Pq  •  <13) 

In  Secs.  2  and  3  V^^Cq)  is  compared  with  the  perturbation  V  .  In 
Sec.  4  approximations  to  veff(<l)  l°r  helium  atom  are  computed. 

It  is  our  hope  that  the  function  V^ff(q)  will  lend  itself  to  empiricism. 


2.  Truncated  Reaction  Operators  and  the  Approximate  Solution  of  the 
SchrPdinger  Equation. 

In  this  section  we  limit  consideration  to  trial  functions  which 
may  be  expanded  in  terms  of  a  truncated  basis  set,,  (nj  s  of  n 
eigenfunctions  of  an  Hq  <>  For  this  choice  of  trial  function  t  is 
equivalent  to  a  truncated  reaction  operator^  t(n)  ^  which  involves 
only  the  set  {%}  .  (The  use  of  an  arbitrary  truncated  basis  set^ 

jcOjJ  ,  is  discussed  briefly  in  Sec.  3;  but  the  treatment  is  more 
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complicated.)  The  matrix  elements  of  t(n)  are  obtained  by  an  itera¬ 
tive  procedure^  after  which  approximations  to  the  energy r  E  (n)  ^  and 
the  normalized  wave  function^  'Wn  (n)jcan.  be  easily  calculated.  Finally 
an  expression  for  the  approximate  “ef fective”  perturbation,  veff^>n^  > 
is  obtained. 

For  the  class  of  trial  functions  which  are  linear  expansions  in 
the  truncated  set  m  ,  the  operator  P  is  equivalent  to 


p<">  -  £ <?v 


(14) 


k=l 


with  0  4*  P(n)  “  1  »  Using  Eq0  (14)  we  may  write 


To(n)  -  P(n)/(E  (n)  -  Ho> 


(15) 


where  the  inverse  operator  in  Sq.  (15)  is  well  defined  assuming  E  (n) 

(o)  q 

does  not  equal  any  of  the  eigenvalues  .  Using  Eqs.  (1),  (3),  (6), 

(14)  and  (15)  we  obtain 


n 

t(n)  -  V  +  V^j^Xg 
k=>l 


k> 


(16) 


.  ,  (E  (n)  -  E,(o)) 
k^q  k  ' 


and 


y»>  ■  K$°>+  <% It<»)l90,> 


(17) 


The  approximate  wave  operator.,  vf(n)  }  may  be  obtained  from  Eqs. 
(9)  and  (16); 


W(n) 


(I/V) t(n) 


(18) 


Using  Eqs„  (8)  and  (18) s  the  normalized  approximate  perturbed  wave 
function  is  given  by 


n 


k»l 


(19) 
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where 


Ckq(n)  =  Ntkq(n)/(Eq(n)  “  Ek  > 


n 

■  a  t 


2  ■)-% 


>n 


Thus  both  Eq(n)  and  Ajy  (n)  are  easily  obtained  once  the  tfcq(n) 
are  computed. 

The  matrix  representation  of  t(n)  }  using  Eq.  (16) ,  is 
n 

-  IvWn)/  (Eq0>  +  V">-*k°>>  > 

k=1  (22) 

The  system  of  equations  given  by  Eq.  (22)  may  be  solved  iteratively.  A 
guess  is  made  for  tqq(n)  •  Then  the  s'ystem  of  equations 


t  (n) 
sq 


11 

■r 


T.kV0,<*?>  +  Vn)  •Ek°))  ' 


s=l^ o o o 9 q-l^q+1^ 0 • • ^n  , 

is  solved  for  the  tfc  (n)*s  with  k*q  .  These  tkq(n)’s  are  then 
used  in  the  expression 


n 

,«  -  E  Vi.q(n)/<E$0)  +  V0  '  Efc0>) 


to  calculate  a  new  value  for  t  (n)  •  This  procedure  is  repeated  until 
self  consistency  is  obtained.  From  the  final  values  for  the  tkq(n^  3 
the  approximate  perturbed  energy  and  normalized  wave  function  are  given 
by  Eqs.  (17)  and  (19)  respectively5.  Similar  iterative  methods  are 


5  Solutions  may  also  be  obtained  for  the  reaction  operator  matrix 

elements  tk«(n)  ,  with  A  /q  ,  using  Eq.  (16).  These  matrix  elements,, 
however,  do  not  have  any  obvious  application  in  the  formulation  for 
non-degenerate  states.  Here  it  should  be  remembered  that  the  t-opera- 
tor  is  defined  with  respect  to  the  state  under  consideration,  q  . 
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available  for  solving  the  secular  equation  of  the  linear  variational 
method^  where  the  variational  coefficients  are  treated  in  much  the  same 
way  as  the  t^^(n)  (see  Sec.  3  and  for  example  ref.  6)  . 

The  connection  between  the  approximate  "effective”  perturbation., 
Veff(^jti)^  and  V  is  obtained  using  Eqs.  (11)  and  (16); 

n 

W’>">  ■  ■  w9='q)Itk,fk/(V°)  • 

k'1  (25) 

This  function^  ■>  equivalent  to  the  operator  t(n)  in 

the  computation  of  the  energy  and  wave  function  of  the  perturbed  system 
(see  Sec.  1);  indeed  t^^(n)  *  (Veff «  An  approximate  "effective" 
perturbation  is  computed  for  a  specific  example^  the  He  atom_,  in  Sec. 

4.  In  Sec.  3  an  expression  for  Veff(q..,n)  is  obtained  in  terms  of  an 
arbitrary  basis  set  -f 60^- 

3.  The  Reaction  Operator  and  the  Variational  Method 


The  reaction  operator  formulation  is  closely  connected  to  the  linear 

variational  method  for  solving  the  Schrtfdinger  equation.  This  property ^ 

3 

is  of  course^  implicit  in  the  derivation  of  the  reaction  operator.  . 

From  Eq.  (22)  one  obtains 
m 

X  [v  -  <V>  •  Ek0>)  tk,(”>/(Vn)  •  Eko))  -  °  * 

k=1  (26) 

■  ‘  -  s=l j 2  y  ° • *  > n 

The  matrix  element  of  the  perturbed  Hamiltonian^  H  ,,  is  given  by 


sk 


=  (H  + 
o 


v) 


sk 


= 


k  °ks 


+  V 


sk 


(27) 


Hence  Eq.  (26)  becomes 


P.  0.  LCwdin^  J..  Molecular  Spectrosc.  10 ,  12  (1963). 
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I[Hsk-  Vn)Sjtk,<")/(Vn)  '^o))  ■  0  >  <28> 


This  set  of  equations  has  a  non~trivial  solution  only  if  the  deter¬ 
minant  al  equation 


sk> 


=  0 


(29) 


where 


Dsk  -  (\k'V”>^KVn)  ■Ik°>>'1 


(30) 


is  satisfied,,  The  factor  (E^(n)  -  E^0^)  ^  ;  common  to  each  member  of 
the  k-th  column,  may  be  divided  out  of  the  determinant  given  by  Eq0 
(29)*  Thus  we  obtain 


■J-1 


(E  (n)  -  E 

ci 


'  sk 


E0(n)  g 


ks  1 


=  0 


(31) 


Since  the  product  factor  is  not  zero,  Eq«  (31)  reduces  to  the  usual 
secular  equation^ 


-  En<n> 


-  0 


(32) 


The  relationship  between  the  linear  variational  coefficients, 
C^^(n)  ,  and  the  matrix  elements  of  the  reaction  operator,  t^^(n)  , 
is  given  by  Eq„  (20),  This  formula  is  valid  only  for  a  basis  sat 
composed  of  orthonormal  eigenfunctions  of  Hq  . 

Arbitrary  Truncated  Basis  Set.  Consider  an  orthonormal  truncated 
basis  set,  9  t*ie  sing^e  restriction  that  OJ ^  -  (p ^  . 

Using  this  basis  set  we  may  write 


J.  0.  Hirschf elder,  C0  F„  Curtiss  and  R,  B,  Bird,  "Molecular  Theory 
of  Gases  and  Liquids”  (John  Wiley  and  Sons,  Inc, /New  York,  1954) 
p.  63, 
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n 

r*» 


t(n)  -  V  +  v  )  l(jJk>  <OJk|  t(n) 


(33) 


k-1  (E(n)  -  H 
,  ,  q  o 

k#q 


Here*  however*  we  may  not  replace  by  as  was  done  previously 

when  the  set  was  usedo  On  this  account^  there  is  no  simple 

theory  such  as  that  developed  in  Sec0  20  On  the  other  hand  it  is  easy 

to  obtain  the  relationship  between  the  C  (n)  and  the  matrix  elements 

icq 

of  the  reaction  operator. 

Let  us  expand  the  normalized  wave  function*  *\J/^(n)  j  terms 
of  the  basis  set  * 

n  n 

•x^(n)  -  »*q(n>  -  .1^  -  £**<*“*  <34> 

k=l  k»l 

where  N  is  a  normalization  constant.  Using  Eqs„  (8),  (18)  and  (34) 
we  obtain 


(35) 


(36) 


n 

t(n  )<pq  =  Vfq(n)  -  V^(Ckq(n)/N)6t4  . 

k=l 

Thus  the  general  expression  for  the  t  (n)  in  terms  of  the  G^^(n) 
becomes 

n 

tsq(n)  -  ^(Ckq(n)/N)Vsk  ,  s-1,2,.  • .,n  . 
k=a 

Using  Eq.  (35)  one  may  obtain  an  expression  for  n)  ; 

n 

Ve££(,,„)  -  (V/<jDq)^<  Clq(n)/K)04  . 

k«l 

Equations  (36)  and  (37)  contain  Eqs.  (22)  and  (25)  as  special  cases. 

g 

4.  An  Example,  the  Helium  Atom  ;  Discussion. 

We  now  consider  a  simple  example  to  illustrate  some  of  the  ideas 


(37) 


A  more  detailed  version  of  this  work  is  available  in  the  technical 
report*  WIS-TCI-10*  Theoretical  Chemistry  Institute*  University  of 
Wisconsin*  25  February  1963. 
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developed  in  the  preceding  part  of  this  paper.  The  truncated  reaction 

2  2  1 

operator  matrix  elements  are  obtained  for  the  (Is)  5 (Is) (2s) 9 (2s)  S 
states  of  helium.  Simple  basis  sets  are  used,,  which  consist  of  one  or 
more  configurational  interaction  wave  functions  <>  constructed  from 
products  of  uniformly  scaled  hydrogen  Is  and  2s  orbitals.  The  approxi¬ 
mate  "effective”  perturbation,,  V^fCq^n)  9  is  computed  for  specific 
examples.  Finally,,  the  results  obtained  for  helium  and  the  more  general 
results  of  the  previous  parts  of  this  paper  are  discussed. 

The  two  electron  problem  is  very  advantageous  for  this  investiga¬ 
tion  because  of  simplifications  resulting  from  the  factorization  of 

electron  spin.  This  factorization  of  spin  permits  the  use  of  the 

9 

following  basis  functions  ; 


where 


II 

ls(l) 

ls(2) 

Vl  - 

(2)"4 

[ls(l)  2s(2)  +  Is (2)  2s(l)J 

£>3  = 

28(1) 

2s(2) 

ls(i) 

=  GO' 

(z)3/r2exp(-Zri) 

2s(i) 

=  (4)“ 

1  (2jt) (Z)3/2(2  -  Zr.)exp(-Zri/2) 

3  and 

Z  is  a 

i  scaling  parameter  which  may  be  energy 

(38) 


(39) 


The  Hamiltonian^  H  ^  for  the  helium  atom  may  be  written  in  the 
form  of  Eq0  (2)  with 


Ho  -  »(%)V22  -  z  [d/rp  +  (l/r2)] 


(40) 


9 

All  equations  and  data  given  in  this  paper  are  in  atomic  units^ 
where  the  unit  of  length  is  a  9  the  Bohr  radius^  and  the  unit  of 
energy  is  e^/a0  • 
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V  -  (Z  -  2)  [(l/rp  +  (l/r2)]  +  l/i^ 


The  eigenvalues  Ev 


1>,2j)3<9  are  easily  computed  by  applying  H 


to  the  functions  given  by  Eqs„  (38 )  •  The  matrix,  elements  of  the  pertur¬ 
bation,  V  ,  can  be  computed  using  the  tables  of  Roofchan^  and  Barnett 
and  Coulson^o 

To  obtain  the  truncated  reaction  operator  matrix  elements  for  the 

2  2  1 

(Is)  , (Is)  (2s)  , (2s)  S  states  of  helium  we  set  q  ®*  1,2,  and  3, 

respectively,  in  Eqs*  .(22),  (23)  and  (24) ,  and  use  the  iterative  method 

12 

of  solution  outlined  in  Sec0  2  «  In  the  Appendix  we  discuss  briefly 

the  selective  nature  of  the  iterative  method  of  solving  for  the  reaction 
operator  matrix  elements,,  It  is  shown  that  the  procedure  will  converge 
(for  the  calculations  considered  in  this  paper)  to  the  state  under  con¬ 
sideration  regardless  of  the  initial  choice  of  t  (n)0  In  Table  I  (page  11)  the 
t^(n)  are  listed  at  several  stages  of  approximation  '  «  Each  approxi¬ 
mation  is  computed  using  one  or  more  of  the  configurations  (Is) (Is) (2s) , 
2 

(2s)  in  the  basis  set„  For  the  unsealed  computations  with  Z  «  2  , 

the  t^(n)  decrease  by  approximately  an  order  of  magnitude  as  k 

takes  on  the  values  1, 2,  and  3  respectively «  Similar  trends  for  the 

t^^(n)^  with  q  /  1^  are  suggested  for  k  >  q  (see  Table  I) „ 

However  with  the  limited  data  given  here  it  is  impossible  to  establish 

definite  trends  for  these  matrix  elements,,  The  fc^(n)  9.  for  Z  -  2  f 

seems  to  be  a  relatively  slowly  varying  function  of  n  „  Similar  but 

less  pronounced  trends  in  the  t,  (n)  are  seen  in  the  energy  optimized 

kq 

computations  0 


C.  C„  J0  Roothan,  J.  Chem0  Phys.  19,,  1445  (1951) . 

M.  Po  Barnett  and  0o  A0  Goulson,  Phil.  Trans 0  (Series  A)  243  s  221 
(1951). 

The  iterative  method  was  carried  out  numerical] y  using  a  Control 
Data  Corporation  1604  computer,,  In  the  energy  minimization  part  of 
the  program  a  simple  search  method  was  employed,,  All  data  has  been 
rounded  off„ 


The  values  of  tk^(n)  ,  with  JL  /*  q  ,  are  not  given  in  this  paper 
(see  refo  5),  but  are  available  in  Tables  IV,  V,  and  VI  of  ref„  8„ 
The  physical  significance  and  use  of  these  matrix  elements  is  at  the 
present  time  unclear,, 


Having  computed  the  truncated  reaction  operator  matrix  elements 

one  may  compute  approximations  to  the  normalized  perturbed  wave 

2  2  1 

functions  for  the  (Is)  (Is)  (2s)  9  (2s) '  S  states  of  helium  using  Eqs„. 

( 19)  ^  ( 20)  j,  and  (21);  these  are  given  in  Table  II  (page  13).  For  both 
the  unsealed  and  energy  optimized  calculations  the  trends  in  the  numer¬ 
ical  values  of  the  coefficients  of  the  normalized  perturbed  wave 

are  more  clearly  defined  than  those  for  the  t^^(n) 
(compare  tables  I  and  XI).  The  G^^(n)  show  a  much  smaller  change  in 


function^,  (n)  i 


magnitude  upon  scaling  than  the  t^^(n) 

C7  (n)  act  as  a  "normalized"  tT  (n)  . 
kq  ,(o)kq 


o  Thus  we  conclude  that  the 
This  effect  is  due  to  the 


energy  difference^  E^(n) 


appearing  in  the  definition  of  the 


C,  (n)  given  by  Eq0  (20).  This  type  of  trend  in  the  coefficients  of 

kq 

a  linear  variational  wave  function  is  aj->parent  in  the  work  of  others 
--  for  example,,  the  Shull  and  LI5wdin  investigation  of  the  ground 


state  of  helium' 


14 


With  numerical  experience  for  atomic  and  molecular 


problems  one  may  eventually  be  able  to  empiricize  the  t.  (n)  orA 

kq 

equivalently  <>  the  C,  (n)  for  helium  and  other  systems 0  Good  estima~ 
kq 

tions  of  the  energy  could  then  be  obtained  from  Eqs0  (17)  and  (36) 0 

One  can  easily  compute  approximations  to  the  energy  of  the 
2  2  1 

(Is)  j (Is) (2s) s (2s)  S  states  of  helium  using  Sq0  (17)  and  the 
t  (n)  obtained  by  the  iterative  procedure  of  Sec.  2.  These  energies , 
given  in  Table  III  (page  14) *  are  not  particularly  good  because  of  the 
nature  of  the  truncated  basis  set  used. 

The  approximate  "effective"  perturbation 9  (q(,n5Z)  s  can  be 

computed  using  the  data  of  Table  I  (or  Table  II)  and  Kq.  (25)  (or  Eq0 
(37)).  As  an  example  we  give  Vpff(1^3^2)  and  )  where 

Z  t^ie  ener§y  optimized  Z  .  Defining 


g(i«,a)  =*  (1  -  ar.)exp(ar^)  «,  1=1  or  2 

V(Z)  «  (Z  -  2)  Rl/rp  +  (l/r2)]  +  1/r 


12 


(42) 

(43) 
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L3 


I 


(Is) (2s)  -0.302599  0.952848  0.022684  -0.1226  0.9923  -0.0570 

(2s)2  -0.016178  -0.028931  0.999451  -0.0239  0.0999  0.9983 


I 


a 

§ 

H 


functions 

CM 

CO 

W 

1 

-0.72184 

00 

rM 

rM 

CN 

r^ 

o 

n 

CM 

O 

<t 

rM 

n- 

o 

D 

cn 

in 

cn 

rH 

n> 

o 

0 

05 

■■■ 

CS 

£ 

H 

H 

H 

a 

d 

■ 

d 

d 

£ 

*rl 

a) 

♦rl 

IT) 

M0 

CM 

CM 

r— 1 

d 

4-J 

CO 

cn 

O 

cn 

0) 

d 

Pt 

H 

in 

id 

rC 

rH 

O 

CM 

H 

m 

H 

lO 

m 

O 

w 

o 

H 

■ 

o 

o 

o 

4H 

d 

P 

• 

H 

■ 

O 

0 

0 

CM 

■ 

■ 

CM 

CM 

CM 

MM 

1 

H 

8 

0 

0 

05 

03 

H 

CL) 

£ 

■■ 

[ 

1  «S 

■P 

o 

60 

HHMCi 

d 

•rM 

P 

m&m 

■ 

■ 

4J 

4-J 

CD 

■ 

05 

CCS 

d 

p 

w 

cO 

O 

m 

lA- 

w 

d 

\0 

■ 

cn 

CM 

<t 

-H 

60 

r> 

00 

<± 

m 

*r-l 

r— 1 

<t  ! 

■ 

■ 

<f 

m 

in 

N 

m 

w 

0° 

■ 

CO 

oo 

■ 

00 

AN 

d 

• 

■ 

© 

O 

o 

05 

o 

CM 

CM 

CM 

CM 

CM 

a 

1 

■ 

H 

U 

0 

O' 

A 

o 

1 

B 

|  H  \  B 

05 

p 

CM 

4-) 

NA 

CD 

■ 

Pj 

■ 

■ 

w 

S 

■ 

a\ 

i—i 

s 

oo 

vO 

o 

w 

03 

'A 

■ 

CM 

CM 

MO 

A 

h- h 

Cn 

oo 

00 

r> 

M 

> ^ 

cn 

IH 

cn 

cn 

cn 

cn 

H 

W 

VO 

M0 

CD 

vO 

03 

T— J 

0 

a 

d 

rM 

CCS 

o 

o 

o 

O 

NA 

■H 

l 

8 

8 

1 

d 

*  •  AiS 

h- h 

cd 

■ 

CD 

p 

CM 

MH 

■m 

MKXH 

rd 

■ 

■  H 

4J 

MM 

II 

■ 

O 

■  SB' 

■ 

h 

P 

CM 

rM 

.■■•■  H 

cn 

CM 

CO 

O 

d 

in 

H 

r— < 

M0 

MM 

d 

cn 

cn 

r-- 

cn 

•H 

P 

■ 

iO 

■ 

■ 

in 

MD 

M0 

A 

0 

CM 

■ 

cn 

'PA'. ."•'  •’•  ■■: 

■ 

in 

cn 

MM 

w 

o 

■ 

B 

cn 

o 

mEM 

As 

d 

■ 

© 

H  H 

• 

• 

a 

d 

CN 

H  H 

rM 

CM 

V 

CCS 

60 

J 

| 

0 

a4 

P 

||| 

■ 

■ 

W 

p 

CD 

m 

.  i 

CD 

d 

w 

Hi 

H 

B 

■l 

A 

£ 

■. 

■ 

■ 

■ 

05 

d 

■ 

a) 

d 

o 

cn 

00 

^■1 

VO 

*H 

o 

1 

■ 

cn 

cn 

■ 

cn 

60 

CD 

o 

cn 

<t- 

cn 

P 

4d 

■ 

o 

■ 

o 

o 

■ 

O 

CD 

4J 

lO 

H 

in 

cn 

^■1 

cn 

d 

■ABi 

r- 

■ 

r-- 

CO 

■ 

co 

CD 

d 

O 

H 

o 

o 

■ 

o 

d 

■ 

CM 

H 

H 

CM 

CM 

HI 

CM 

u 

ccJ 

1 

■ 

8 

0 

HI 

0 

0 

K 

■ 

HI 

4-J 

CCS 

A 

wmm\ 

■1 

IHI 

p 

CM 

(D 

P< 

CM 

o 

A 

An 

P 

05 

d 

o 

CM 

o 

4J 

•H 

O 

A 

pi 

CCS 

05 

A 

CM 

AN 

o 

MM 

a 

f/3 

/~s 

05 

d 

o 

CM 

03 

CM 

CD 

CD 

CM 

v_x 

CM 

'w' 

u 

i — 1 

4-> 

/A 

A 

An 

n 

CCS 

C/5 

r/3 

A 

05 

d 

o 

P 

d 

AN 

CM 

rM 

AN  ‘ 

rM 

CD 

05 

d 

CD 

03 

'w' 

05 

N-*' 

4-J 

60 

03 

CM 

A 

A 

CM 

A 

rt 

CD 

•H 

d 

CM 

w 

CM 

CM 

CM 

CM 

o 

rC 

MM 

/A 

/A 

/  \ 

/A 

A 

AN 

An 

d 

4-J 

d 

03 

03 

05 

03 

03 

05 

05 

d 

0 

rM 

rM 

CM 

rM 

rM 

rM 

rH 

p 

MM 

a 

na 

n_/ 

v _ ' 

_ _ ' 

N— ✓ 

H 

O 

14 


i 


15 


we  obtain 


Veff(l,3,2)/V(2)  =  {l-(0. 079382)  [g(l,l)  +  g(2,lj  -(0.002122)g(I,l)g(2,l)} 


fl-(0. 03090)  fg(ls0o  8792)  +  g(2.,0. 8792)j 
^ff^1*3*1*7583^^1-7583)  =  | -(0. 00302)  g(l, 0o 8792)  g(2>0.  8792) 


The  simple  form  of  the  approximations  to  V  ^  for  the  helium  atom 
given  above  is  due  to  the  limited  truncated  basis  set  used.  The  use  of 
better  basis  sets  would  permit  accurate  calculations  of  V  ^  •  After 
sufficient  numerical  and  theoretical  experience  we  hope  it  will  be 
possible  to  accurately  empiricize  V  ^  for  many  atomic  and  molecular 
problems. 

One  of  the  most  useful  aspects  of  the  reaction  operator  formula¬ 
tion  is  that  the  operator^  t  9  takes  into  account  all  interactions  of 
the  system  above  the  unperturbed  state.  Thus  t  and  its  matrix  elements 
(and  of  course  V  ^  and  Its  matrix  elements)  must  give  valuable  infor¬ 
mation  concerning  the  physical  properties  of  the  perturbed  system. 
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Appendix.  Convergence  Properties  of  the  Iterative  Method  of  Solving 
for  the  Reaction  Operator  Matrix  Elements . 

In  the  computations  outlined  in  Sec0  40,  interesting  convergence 
properties  of  the  iterative  method  for  obtaining  the  reaction  operator 
matrix  elements  were  apparent.  The  method  converged  to  the  t  (n) 
corresponding  to  the  state  under  consideration  even  if  the  initial  guess 
for  tqq(n)  corresponded  to  the  solution  for  another  perturbed  state. 

In  the  following  discussion  we  make  use  of  a  theorem^  which  states: 

If  the  equation,,  x  =  f(x)  subject  to  the  condition  f’(x)  ^  0  ,  is 
solved  by  the  iterative  procedure,  x^+ ^  “  fCx^.)  *  then  the  method 
will  converge  to  a  particular  root  A  *  if  lf,(  A)[<  i  •  if 
I  f *( A) [ >  1  )  the  iterative  procedure  will  diverge  from  the  root  X  . 

Q 

Using  Eq„  (22)  with  q  =*  1;  Z  ~  2;  s  =  1,2,  or  3^  one  can  show 
that  for  the  ground  state 


tu<3)  -  v11  +  vj2  («) 

+  VU  [Yu  -  +  2V12V23V31  |Vl2  Y13 ' 4] 

where 

Tpj  -  e5°}  -  Hjj  +  tn(3)  4  j=233  .  (46) 


If  is  added  to  both  sides  of  Eq.  (45)  we  obtain  the  Feenberg 

1  15 

perturbation  series  for  a  basis  set  comprised  of  three  real  ortho¬ 
normal  eigenfunctions  of  Hq  .  Equation  (45)  has  the  form  t^  =  f ( tz ^ 
Using  Eq.  (45)  one  obtains: 


ff(tn(3))  = 


9  9  9  9  9  2 

V12V23  +  V13V23  +  V12^13 

( - 

CO 

CM 

+  V13^'l2+ 2V13V32V21^12  +  Tl3^ 

CO 
CM  CM 
> 

1 

.  - 

L  J 

(47) 
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17 


The  roots  of  Eq«  (45) }  t£^(3)  }  which  correspond  to  the  j-th 

perturbed  state  can  be  calculated  by  the  following  equation: 


=  t.O) 

-  ,  j=l,2,3  • 

(48) 

The  values  for  |f*(t£|  (3))| 

j  j  =  1^2,3 j  are  given  below: 

if*{tii)<3))i  ■ 

0.1011 

IfUt^o))!  = 

9.960 

(49) 

|f'(t<3)(3))[  = 

1.71  x  103 

Equations  (49)  explain  the  selective  nature  of  the  iterative  procedure 
for  the  ground  state  calculations  with  Z  =  2  n  «  3.  Similar 

g 

explanations  can  be  given  for  the  other  computations  of  Sec.  4. 


